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Abstract 

We demonstrate a standing-wave light pulse sequence that places atoms into a superposition of 
wavepackets with precisely controlled displacements that remain constant for times as long as 1 
s. The separated wavepackets are subsequently recombined resulting in atom interference patterns 
that probe energy differences of ~ 10~ 34 Joule, and can provide acceleration measurements that 
are insensitive to platform vibrations. 
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I. INTRODUCTION 



Atom interferometry employs the interference of atomic de Broglie waves for precision 
measurements jj. In practice two effects limit the ultimate sensitivity of devices where 
the interfering atomic wavepackets are allowed to propagate in free space: the effect of 
external gravitational fields upon the atomic trajectories, and transverse expansion of the 
atom cloud. By accounting for gravity, atomic fountains can increase the interrogation 
time during which the interferometry phase shifts accumulate [2|; alternatively one can use 
magnetic dipole forces to balance the force of gravity Magnetic waveguides 0, 0] can 
trap atoms for times longer than a second, suggesting the possibility of measuring energy 
differences between interfering wavepackets with an uncertainty < h/(l s)~ 10~ 34 Joule; 
however, this remarkable precision cannot be obtained if the decoherence time of the atoms 
is much shorter than the trap lifetime. Early atom interferometry experiments using atoms 



confined in magnetic waveguides showed that the external state coherence o: 
decayed quite quickly, limiting interferometric measurements to times <10 ms 



the atoms 
7]. More 



recent experiments using Bose condensates [8] have shown that the external state coherence 
can be preserved for approximately 200 ms, where the decoherence is dominated by atom- 
atom interactions. Interferometry experiments using either condensed atoms in a weak trap, 
or using non-condensate atoms in a waveguide with precise angular alignments have been 



shown to have phase-stable interrogation time o: 



12j. 



50 ms, where the dephasing is induced 



by inhomogeneities in the confining potential. {9- 

This work demonstrates a new atom interferometer configuration that measures the dif- 
ferential phase shift of spatially-displaced wavepacket pairs. We demonstrate phase-stable 
interferometry operations with up to one second interrogation time by applying the tech- 
nique to atoms in a straight magnetic guide. We show that the matter-wave dephasing rate 
scales linearly with the wavepacket displacement, suggesting that dephasing in our inter- 
ferometer is primarily caused by a weak longitudinal confinement of the atoms. We also 
demonstrate that the phase readout of the interferometer is less sensitive to vibration than 
conventional interferometery schemes, which should enable precision measurements even in 
noisy environments such as moving platforms. 



2 



d 2 - 











MM 









T-2T. 



FIG. 1: (Color online) Recoil diagram of the 4-pulse scheme. The light shift potential of the 
standing wave is represented by the array of circles. Straight lines show the centers of the diffracted 
wavepackets. A representative pair of interfering paths is marked with thick lines. The matter- 
wave interference at around t = T is probed with E a . The backscattered light, e.g., the "grating 
echo" E b (r), is mixed with a weak local field from E b mode for phase retrieval. 



II. THE 4-PULSE GRATING ECHO SCHEME 



Typical Talbot-Lau matter-wave interferometery [13Hl6j employs a 3-grating diffraction 
scheme. In the most common time-domain setup, an atomic wavepacket is diffracted by 
a periodic potential, applied briefly at time t = 0, into a collection of wavepackets that 
depart from each other at multiples of the velocity vq = hQ/m, where Q is the potential's 
wavevector and m is the atomic mass. The potential is pulsed on again at t = T/2, and the 
different velocity classes created by the first pulse, which have now moved away from each 
other, are each again diffracted into multiple orders. After the second pulse there will be 
pairs of wavepackets whose relative velocity has been reversed from before; these will move 
toward each other, then overlap and interfere near time t = T. Those with relative velocity 
nvq will generate density fringes with wavevector ±nQ. 



As in earlier experiments 
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171 ] we use an off-resonant optical standing wave (SW) 
to create the pulsed periodic potentials, and observe the resulting atomic density fringes 
by measuring the Bragg scattering of an optical probe. In our experiment we observe the 
lowest order fringes, n = ±1, where the Bragg condition corresponds approximately to 
backscattering of one of the beams that forms the standing wave. 

The interferometric technique presented here employs a 4-pulse scheme (Fig. [1]), where 
the additional pulse is used to halt the relative motion of the interfering wavepacket pairs. 



3 



After the first pulse is applied the situation is identical to that in the 3-pulse case, with the 
original wavepacket split into different diffraction orders corresponding to velocities ±?ivq. 
Here we quickly apply a second pulse after a short time t — T s . The pairs of wavepackets 
that we are interested in, those that will eventually interfere, are those that now have zero 
momentum difference; these pairs have the same velocity but are displaced in space, having 
moved apart by a distance (d = vqT s ) during the interval between pulses. After waiting 
for a time T — 2T$ the coherence between the wavepackets is measured by allowing them to 
interfere; a third pulse diffracts the wavepackets at t = T — T s and the resulting interference 
fringe is probed around t = T. This 4-pulse sequence can be imagined as a 3-pulse sequence 
of length 2T S that is "paused" between times T = T$ and T — Ts; though during this time 
the relative phase of each separated wavepacket pairs will continue to be sensitive to external 
fields. 

The existing theory of Talbot-Lau interferometry can be straightforwardly extended to 
this scheme. We consider the SW field formed by the traveling light fields E a and E b with 
associated k-vectors k a and k fe , with Q = k a — k 6 = Qe x . Experimentally, we measure the 
backscattering of light from k a into k b ; this is characterized by an electric field component 
E b {r) that can be expressed in terms of the atomic density operator: 

E\t) = -igE p Tr[p(T)e iQ ^}. (f) 

Here E p gives the amplitude of probe light from the E a mode, g is a constant that depends 
on the atomic polarizability and the number of atoms participating in the interaction, x and 
p are the position and momentum operators for atomic motion along e x , and p(T) is the 
single atom density matrix at time T. Consider an atomic sample with a rms velocity u 
and a thermal de Broglie wavelength l c = h/(mu), in Eq. (jTJ) we use r ~ I c /vq to specify 
a coherence- length-dependent time window (T — I c /vq,T + I c /vq) during which the atomic 
wavepackets overlap so that the interference fringe contrast is non-zero. Experimentally, 
the amplitude of E b (r) is averaged during this time window to extract the magnitude of 
the interference fringe; this will be referred to as the amplitude of the "grating echo", 



E 9 (T S ,T):=< E\r)> 



(18|. 



For a SW pulse in the Raman-Nath regime, where the atomic motion can be neglected 
for its duration, the n th order matter- wave diffraction is weighted by the amplitude i n J n (8), 
with J n the n th order Bessel function and 9 the time-integrated light shift or pulse area. We 



4 



define {Ti, T 2 , T 3 , T 4 } = {0,T s ,T — T S ,T}, and specify the position of standing wave nodes 
at Ti with the SW phase (fi. The interaction of the first three SW pulses in the 4-pulse 
interferometer can be effectively described by 

Since the standing wave phases ifi involve simple algebra, we will ignore them during the 
following discussion, and reintroduce them when they become relevant. 

In order to account for imperfections in the guide we consider the ID motion of atoms 
along e x during the interferometry sequence to be governed by H = H Q + Vsw(t) and further 
Hq = p 2 /2m + V(x) where V(x) is a general ID potential. We introduce the time-dependent 
position and momentum operators x(t) = e iHot xe~ °* and p(t) = e lH o t p e -' lH oi _ p or r p s i 
be sufficiently short, the atomic motion can be treated as free during < t < T — T s and 
T — T s < t < T . For a thermal atomic sample with Ql c « 1, we find at the leading order 
of Ql c , the interferometer output is related to the initial atomic density matrix p(0) by 

E g {T s ,T) = -igEpJ-x (20 3 sin u Q T s )x 
E c m Tr[p(5x m ,5p m ,0)e ± n d ^e-^% ^ 

mi ,ni2 

Here c m = J m M Jm 1+ i(9i) JmM J m2 -i(^)e^, <\> m = (2mi + l)Ql e + (2m 2 - l)u Q T a , 
5x m = ttiiVqTs, Sp m = (mi + fn^jQ- 111 Eq. ^ uq is the two-photon recoil frequency of 
atoms and we have d = vqT s . 

The second line of Eq. ()3]) composes a weighted sum of matter-wave correlation func- 
tions. The initial conditions of matter-wave states are specified by a density matrix 
p(5x, 5p, 0) = e l ( 5x P+ s P x ) p e -t{Sxp+Spx) describes an atomic ensemble that is identical 

to p(0), but with mean position and momentum shifted by Sx and Sp respectively. The 
correlation function gives the average overlap of wavepacket pairs propagating under an ex- 
ternal potential displaced by cfe x , with one example sketched with the thick lines in Fig. [TJ 
The correlation functions are in direct analogy to the neutron scattering correlation func- 
tion discussed in ref. where momentum displacements were considered. Notice that 
if the uniform atomic sample has a spatial extension L » 5x and with thermal velocity 
u » Sp/m, the original and shifted density matrix are approximately the same, and the 
correlation functions are approximately independent of 5x m or 5p m . We can thus use a sum 
rule of Bessel functions to simplify the second line of Eq. (J3J) giving: 
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E g (T s ,T) « -igJl^Oz smoo Q T s )x 

where we have chosen 9 2 = 9 3 and define g = gE p 9iQl c . 
Three features of Eq. ([3]) and Eq. @ are worth noting: 

First, though we have only considered the ID motion of atoms in the external potential 
V(x), the formula is readily applicable to a 3D time-dependent potential V(r,t) as long 
as the external potential contributes negligibly to the differential phase shift of wavepacket 
pairs during < t < T s and T — T s <t < T. 

Second, the reduction from Eq. ([3]) to Eq. (pE} requires that the correlation functions 
be insensitive to momentum transferred by the SW pulses. This is very well satisfied if 
the displacements are much smaller than the position and momentum spreadings of the 
atomic gas itself since p(Sx, dp, 0) ~ p(0). In addition, the approximation is particularly 



well satisfied if the potential is periodic at small wavelengths |20l. \21\. 

Finally, notice that both Eq. ([3]) and Eq. (jl]) can be evaluated semi-classically by replacing 
Tr[p(0)e^ T) e-^ dm } with < e^WTJ-^o)] > c , where v(t) gives the classical velocity of 
atoms along e x and < ... > c gives the classical ensemble average over the atomic initial 
conditions. 

Now consider a weak quadratic potential V(x) = max + ^mufx 2 with an acceleration 
force ma and with uiT << 1 to model the potential variation along the nearly free (axial) 
direction of propagation in the magnetic guide, also assuming an atomic sample with a 
Gaussian spatial distribution along e x given by p(x) = e~ x2 ^ 2Ij2 \ The expected amplitude 
of the grating echo signal is then found to oscillate with T s and to decay as a Gaussian with 
the total interrogation time T: 



E g (T Sj T) w -%~g.J\{2B 2 smu Q T s )x 

(5) 

e 2 » 



1 / md Irn t \2 ■/ tnd rn i \ 



(Here we have reintroduced the standing wave phase in Eq. ([2]), where ipij = <pi — ipj) 



III. EXPERIMENTAL SETUP 



The experimental apparatus is described in detail in Ref. 18|. A straight 2D 
quadruple magnetic field with a transverse gradient of 70 G/cm is generated by four 
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FIG. 2: (Color online) Examples of interferometry signal amplitude C(T S ,T) oscillation vs T s at 
fixed T. Scatter plots are from experiments. Solid lines are calculated according to an extension of 
Eq. © with complex pulse area 6 = 9(1 + 0.025i). (a): 9 = 0.85 and T = 50 ms. (b): 9 = 2.8 and 
T = 30 ms. The data were taken during different period of experiments with different signal/noise. 

200 mmxlOO mmxl.56 mm permalloy foils poled in alternating directions. Approximately 
10 7 laser-cooled 87 Rb atoms in the ground state F=l hyperfme level are loaded into this 
magnetic guide, resulting in a cylindrically-shaped atom sample 1 cm long and 170 //m 
wide. The average transverse oscillation frequency of the atoms in the guide is on the order 
of 80 Hz, estimated by displacement induced oscillations of the atomic sample using absorp- 
tion images. A very weak harmonic potential along the guiding direction is estimated to be 
coi ~ 2tt x 0.08 Hz |l8|. 



The SW fields formed by two counter-propagating laser beams with diameters of 1.6 mm 
are aligned to form a standing wave with k-vector along the magnetic guide direction e x . 
Precise angular adjustment is achieved by tuning the orientation of the magnetic guide 
using two rotation stages to within 2 x 10 -4 radians. The optical fields are detuned 120 
MHz above the F = 1 - F' = 2 D2 transition. We choose the SW pulse with typical pulse 
area of 9 ~ 0.8 — 3.0, and with duration of 300 ns to be deep in the thin-lens regime of the 
25 /zK atomic sample. With this pulse duration, the fraction of atoms contributing to the 
final interference fringe is typically limited by SW diffraction efficiency to about ten percent. 
We probe the A/2 atomic density grating at around time t = T by turning on only one of 
the traveling wave beams; the other beam is attenuated and shifted by 6 MHz to serve as 
an optical local oscillator, where the combined intensity is measured using a fiber-coupled 
avalanche photodetector. The beat signal is measured and numerically demodulated using 
the 6 MHz rf reference to recover the grating echo signal E g (T s ,T) = C(T s ,T)e lLp ( Ts ' T \ 
The interferometer signal amplitude C(T S ,T) and phase <p(T s ,T) are measured for different 



7 




C(T S ,T) O) 



<p(T s ,T) (b) 








-350 







100 



200 



300 







100 



200 



300 



T s [ms] 



T s [ms] 



FIG. 3: (Color online) (a): Backscattering amplitude decay. The scatter plots give C(T) at four 
different displacement d. The solid lines are due to Gaussian fit. (b): Interferometry phase shift 
<p(T) vs T at different displacement d. Four of the data traces from the same experiments as those 
in (a) are plotted with thick lines. 

interferometer parameters. 

IV. RESULTS AND DISCUSSIONS 

According to Eq. (jSJ), the pre-factor J\ (26 sin ujqT s ) in the backscattering amplitude is 
an oscillatory function of T s , with the periodicity determined by 2u>q = 27r/33.15 /is -1 . 
The amplitude oscillation is reproduced experimentally; two examples are plotted in Fig. |2l 
where T s is varied from 0.16 ms to 0.26 ms. In Fig. [2(a) a relatively small SW pulse area 
6 ~ 0.85 was chosen so that the Bessel function is approximately linear. Correspondingly, 
we see the oscillation is approximately sinusoidal. In Fig. |2] (b) a strong SW pulse with 
area 9 pa 2.8 was chosen and the Bessel function becomes highly nonlinear. Nevertheless, 
the experimental data still fits the theoretical expectation from Eq. (jSJ) fairly well. The 
values for the of pulse area in the calculation were found in agreement with the SW pulse 
intensity-duration products. Notice the solid lines in Fig. [2] were calculated according to an 
extension of Eq. (jHJ) with a complex SW pulse area including an imaginary part to account 
for the optical pumping effect at the 120 MHz SW detuning [181 ]. 

With fixed T s at the peak values of the amplitude oscillations, we now consider the 
dependence of the interferometer signals on the total interrogation time T. Fig. [3] gives 
examples of the interferometer amplitude decay C(T) and phase shift <p(T) at various d = 
vqT s . From Fig. E(a) we see the amplitude decay is slower for smaller d, while all the C(T) 
fit fairly well to Gaussian decay, in agreement with Eq. (JHJ) derived from a weak harmonic 
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FIG. 4: (Color online) (a): Backscattering amplitude decay rate 7(d) [extracted from a Gaussian 
fit of C(T) data such as those from Fig. [3] (a)] vs wavepacket displacement d. (b): Interferometry 
phase shift rate f(d) = [extracted from a linear fit of <p(T) data from Fig. [3] (b)] vs d. The 

insets give the residuals of the linear fits. 
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FIG. 5: (Color online) Interferometry readouts C(T) and ty?(T) at (f = 0.413 /jm. (a): C(T) vs T; 
(b): 9?(T) vs T, with inset giving the residual after a linear fit. 

confinement model. In Fig. Efb) we see the phase readout is a linear function of interrogation 
time T, also in agreement with Eq. (jSJ). By applying Eq. (j3J) to the observed phase shifts, we 
consistently retrieve an acceleration a = 83.4 mm/s 2 for different d. The acceleration is due 
to a small component of gravity along the standing- wave/magnetic guide direction e x 18], 
as confirmed by varying the tilt angle of the apparatus. 

We extract the amplitude decay rate j(d) by fitting the C(T) decay data with C(T) oc 
e -(7(X)r) 2 _ fj^g dephasing rate 7(d) is plotted vs the displacement d in Fig. 11(a). The 
(^-dependence of 7(d) shows good agreement with a linear fit. According to Eq. fl5]), for 
weak confinement along e x with uji ~ 27r x 0.08 Hz and for L ~ 2.5 mm of our 1 cm 
atomic sample, we expect 7(d) ~ mu fdL/h/y/2 ~ 0.6d/ /im s . This agrees with the 
experimentally measured 7(d) = (2.97 + 0.665<i/ / um) s" 1 according to Fig. H^a). The offset 
of j(d — > 0) = 2.97s -1 is partly due to the escape of atoms from the interaction zone via 
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FIG. 6: (Color online) Interferometry phase readouts (p(T) for a 3-pulse (open triangles, right and 
top axes.) and a 4-pulse (solid squares, left and bottom axes. T s = 0.4018 ms.) configurations. 
Fig. (a) gives the phase readouts. Fig. (b) gives the residual of a quadratic (3-pulse) and a linear 
(4-pulse) fit. 

collisions with the walls of the 4 cm vacuum glass cell, which, if fit to a Gaussian, would 
give j(d — > 0) ~ 1.6 s -1 . The remaining discrepancy is likely due to the inaccuracy of the 
Gaussian fit which is based on the assumption of a weak harmonic perturbation V(x) in 
Eq. (151). For small d and thus a small dephasing rate, local anharmonicity in V(x) might 
become important. Indeed, for long interaction time T the decay exhibits an exponential 
feature, which is clearly seen in Fig.[5](a) where the amplitude decay C{T) with d = 0.418 /zm 
and T s = 35.4 /is is plotted. For such a small wavepacket displacement d, the phase of the 
backscattering signal remains stable for T > 1 s, as shown in Fig. 0(b). 

Last, we consider the effect of phase noise in the SW on the sensitivity of our device, 
induced for example by vibrations of mirrors in the SW path. For T >> T s , the stand- 
ing wave phase variation due to time dependent changes in mirror positions is given by 
<pi,2, which is not correlated with (p^. If we specify the SW phase at time t with (j>(t) 
such that ifi = 0(Tj), the mirror vibration induced interferometer phase noise is given 
by -/Vp jm i rror (T s , T) = < [4>(t) — 4>(t + T s )} 2 > t , which does not depend on T. This is 
different from a 3-pulse atom interferometer with mirror-induced phase noise given by 
A^, mirror (T) = ^ < [4>{t) — 2(f)(t + T/2) + (p{t + T)] 2 > t , where increases in sensitivity due 
to increases in interaction time necessarily also result in increases in phase noise. In con- 
trast, in the four pulse scheme considered here T can be increased to improve the sensitivity, 
while keeping iV v , jmirror (r s ) unaffected. 

This is illustrated in Fig. |6] where we compare the 3-pulse and the 4-pulse interferometer 
phase readouts under the same noisy environmental conditions. A white noise voltage source 
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is filtered to eliminate frequencies below 100 Hz, then amplified and applied to a piezo-driven 
mirror in the SW optical path. As shown in Fig. [61 the mirror vibration randomizes the 
phase of the 3-pulse interferometer for T greater than 5 ms. Under the same conditions, the 
phase of the 4-pulse interferometer is stable for times longer than 150 ms. In this case the 
acceleration sensitivity of the 4-pulse interferometer 5tp/5a ~1 rad/mm/s 2 at T = 150 ms, 
exceeds that for the 3-pulse case of 5ip/5a ~0.4 rad/mm/s 2 at T = 10 ms. The insensitivity 
of the 4-pulse scheme to low-frequency mirror vibrations is a feature of speed-meters, as 
shown in Eqs. @, fl5]) in the semiclassical limit with the phase proportional to the velocity 
during the interrogation time T. 



V. SUMMARY 



We have demonstrated a 4-pulse grating echo interferometer scheme to study the dephas- 
ing effects for atoms confined in a magnetic guide. We find linearly reduced dephasing rate at 
reduced wavepacket displacements, indicating that the matter-wave dephasing is due to very 
weak potential variation along the waveguide in our setup. We have demonstrated phase 
stability for an interferometry sequence with total interrogation time exceeding one second. 
We also show that a four pulse interferometer can provide acceleration measurements with 
very long integration times that are insensitive to apparatus vibrations, though it is impor- 
tant to note that the sensitivity of the interferometer scheme we describe is compromised 
by the small wavepacket separations flo| . 

In the future, such a system could study the quantum stability of wavepackets due to dis- 



placed potentials [lj)| by deliberate 



along the waveguide direction 
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y introducing time dependent variations in the potential 



2lL Instead of measuring the mixed-state correlation 



functions, fidelity-type measurement 19[ can be proceeded with sub-recoil cooled atoms oc- 



cupying a sing 



applied 



22 



23). 



e matter-wave state, where velocity-selective beamsplitting schemes can be 
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